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We critically examine the internal consistency of a set of minimal assumptions entering the the- 
ory of fault-tolerant quantum error correction for Markovian noise. These assumptions are: fast 
gates, a constant supply of fresh and cold ancillas, and a Markovian bath. We point out that these 
assumptions may not be mutually consistent in light of rigorous formulations of the Markovian 
approximation. Namely, Markovian dynamics requires either the singular coupling limit (high tem- 
perature), or the weak coupling limit (weak system-bath interaction). The former is incompatible 
with the assumption of a constant and fresh supply of cold ancillas, while the latter is inconsistent 
with fast gates. We discuss ways to resolve these inconsistencies. As part of our discussion we 
derive, in the weak coupling limit, a new master equation for a system subject to periodic driving. 



I. INTRODUCTION 



The theory of fault-tolerant quantum error correction 
(FT-QEC) is one of the pillars that the field of quan- 
tum information rests on. Starting, with the discovery 
of quantum error correcting codes and the subse- 

quent introduction of fault tolerance , this theory has 
been the subject of many improvements and important 
progress 

iiissiEiEiEiEiEiEiEiEa, 

leading to the well-known error correction threshold con- 
dition. Most recently, work by Steane and Knill 
[l7| (see also Reichardt pUf) has pushed the threshold 
down to values that are claimed to be very close to 
being within experimental reach. A notable feature of 
much of the work on FT-QEC is that the error mod- 
els are phenomenological. By this we mean that the 
underlying models often do not start from a Hamilto- 
nian, microscopic description of the system-bath interac- 
tion, but rather from a higher level, effective description, 
most notably that of Markovian dynamics. E.g., Knill 
writes: "We assume that a gate's error consists of ran- 
dom, independent applications of products of Pauli op- 
erators with probabilities determined by the gate" (our 
emphasis) 17]. This approach is natural given the con- 
siderable difficulty of obtaining error thresholds start- 
ing from a purely Hamiltonian description. Neverthe- 
less, Hamiltonian approaches to decoherence manage- 
ment in a fault-tolerant setting have been pursued, e.g., a 
mixed phenomenological-Hamiltonian treatment of FT- 
QEC [!l E IHEJ, and a fully Hamiltonian study of 
fault tolerance in dynamical decoupling j2^. Also note- 
worthy are recent mixed phenomenological-Hamiltonian 
continuous time treatments of QEC pa. l24l l25 | . 



Here we are concerned with a critical re-evaluation 
of the physical assumptions entering the theory of FT- 
QEC. We scrutinize, in particular, the consistency of 
the assumption of Markovian dynamics within the larger 
framework of FT-QEC. We point out that there may be 



an inherent inconsistency in the theory of Markovian FT- 
QEC, when viewed from the perspective of the validity 
of the Markovian approximation. We begin by briefly re- 
viewing, in Section ^ a set of minimal and standard, 
universally agreed upon assumptions made in Marko- 
vian FT-QEC theory. We then review, in Section IIIII 
the derivation of Markovian master equations, empha- 
sizing the physical assumptions entering the Markovian 
approximation, in particular the requirement of consis- 
tency with thermodynamics. Having delineated the set of 
assumptions entering FT-QEC and the quantum Markov 
approximation, we discuss in Section llVl the internal con- 
sistency of Markovian FT-QEC theory. We point out 
where according to our analysis there is an inconsistency, 
and discuss possible objections. In Section IVII we then 
discuss how one may overcome the inconsistency using 
a variety of alternative approaches, including adiabatic 
quantum computing (QC), holonomic QC, topological 
QC, and recent work on FT-QEC in a non-Markovian 
setting [ljj, |2(J, |2l| ■ We conclude in Section IVTTl 



II. REVIEW OF STANDARD ASSUMPTIONS 
OF FT-QEC 

The following are a set of minimal assumptions made 
in the theory of FT-QEC H S S S S H E IH 13 
El El El E1II3: 



1. Al — Gates can be executed in a time r g such that 



TnUJ 



frequency. 1 



0(tt), where lj is a typical Bohr or Rabi 



1 One might object that slower (even adiabatic) gates could be 
used instead. We analyze this possibility in detail in Section lV Dl 
and show that it does not lead to an improvement. 
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2. A2 — A constant supply of fresh and nearly pure 
ancillas: at every time step we are given a supply 
of many qubits in the state |0), each of which can 
be faulty with some error parameter 77 <C 1. 

3. A3 — Error correlations decay exponentially in 
time and space. 

Some remarks: 

(i) Al is not typically stated explicitly in the FT-QEC 
literature, but can be understood as resulting from the 
definition of a quantum gate, which is a unitary trans- 
formation U = exp(iA); when A — T g H, where H is a 
Hamiltonian generating the gate, Al follows from the 
absence of a free parameter: when r g is scaled up H 
(and hence its eigenvalues) must be scaled down, and 
vice versa. 

(ii) The distinction between Bohr and Rabi frequencies 
in Al is related to the application of constant vs periodic 
control fields, respectively. In the case of a constant con- 
trol field Al can be understood as the condition that sat- 
urates the "Margolus-Levitin theorem" (2(| , which states 
that the time required to transform an initial state 

to an orthogonal state \ip ) using a constant Hamilto- 
nian H is lower-bounded by r m ; n = irh/(2E), where 
E = (ip\H\ip); when is an eigenstate of H we have 
T g > 7r/(2a;), where uj — E/H is the Bohr frequency. See 
also (23 for the adiabatic version of the Margolus-Levitin 
theorem, and [28j for a lower bound on the amount of en- 
ergy needed to carry out an elementary logical operation 
on a quantum computer, with a given accuracy and in 
a given time. In the case of periodic control fields one 
can understand Al as the result of the standard solution 
to the driven two-level atom problem, where the proba- 
bility of a transition between ground and excited state 
is given by (Or/Q^) sin 2 (f2^i/2), where Or, is the Rabi 
frequency and O^ = (Or + S 2 ) 1 ^ 2 , where S is the detun- 
ing. This expression for the transition probability yields 
r g Q' K = 0(n). 

(iii) A2 is shown to be necessary in j6j. A3 is stated 
clearly in pfl (see the discussion in Sections 2.10 and 
10 there). These, and additional assumptions [such as 
constant fault rate (independent of number of qubits) 
and parallelism (to correct errors in all blocks simulta- 
neously)] are explicitly listed, e.g., also in j2^, Section 
II. 

(iv) A3 is usually related to the Markovian assumption, 
however both notions, the space-time correlations of er- 
rors and the Markovian property, need some comments 
and explanations. Using the convolutionlcss formalism 
in the theory of open systems (see, e.g., [3,0] ) it is always 
possible to resolve the total superoperator A(t) as 

n 

A(t) = JjAi0i (!) 

where Ui are ideal unitary superoperators (corresponding 
to quantum logic gates), and Aj are linear maps, not 
necessarily completely positive (CP) or even positive. If 



Aj are CP then we can always realize them by coupling to 
an evironment which is "renewed each time step" . This is 
the "Markovian condition" as formulated in Q (section 
2.10). However, complete positivity is not a necessary 
condition for QEC, which only requires a linear structure 
[3lll32T |. To obtain the Threshold Theorem one needs the 
following bound on the probability 7][Eq. (2.6)]: 

Pr(fault path with k errors) < cr; fc (l - n) v ~ k , (2) 

where rj is the probability of a single error, c is a cer- 
tain constant independent of r), and v is the number of 
error locations in the circuit. This bound implies that 
the fc-qubit errors should scale as ~ n k , i.e., that in the 
decomposition of Aj into fc-qubit superoperators Lj(k) 

\\Lj(k)\\ < cr, k . (3) 

As discussed in |33j (within the Born approximation), 
the condition © can strictly be satisfied only for tem- 
porally exponentially decaying reservoir correlation func- 
tions, while for realistic reservoir models the temporal 
decay is generically powerlike. The decay of reservoir 
correlation functions (i.e., localization in time) translates 
into localization of errors in space due to the finite speed 
of error propagation. On the other hand it is widely be- 
lieved that the Markovian model can be understood as 
arising, to an excellent approximation, from coupling to 
a reservoir which is not only renewed at each time step, 
but whose influence is independent of the actual Hamil- 
tonian dynamics of the open system, and is localized in 
space (independent errors model) |34| . A large part of 
the present paper is devoted to a critical discussion of 
this claim. 

(v) We note that the recent papers on FT-QEC theory 
[13. 120I l2l| relax the (Markovian) assumption A3, but 
do make Al (implicitly) and A2. We comment on these 
papers in Section fVl El 

III. REVIEW OF MARKOVIAN MASTER 
EQUATIONS 

The field of derivations of the quantum Markovian 
master equation (MME) is strewn with pitfalls: it is in 
fact non-trivial to derive the MME in a fully consistent 
manner. There are essentially two types of fully rigorous 
approaches, known as the singular coupling limit (SCL) 
and the weak coupling limit (WCL), both of which we 
consider below. See, e.g., the books [3|j, |3(| for more 
details, as well as the derivation in [37|. 

Consider a system and a reservoir (bath), with self 
Hamiltonians Hg and H R , interacting via the Hamilto- 
nian Hsr = XS ® i?, where S (R) is a Hermitian system 
(reservoir) operator and A is the coupling strength. A 
more general model of the form Hsr = XL ^a^a R a 
can of course also be considered and results in the same 
qualitative conclusions. Thus the total Hamiltonian is 

H =(H Q s +H c (t))®I R + I s ®H R + H SR , (4) 
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where Hc(t) describes control over the quantum device 
(system), and I is the identity operator. 

The SCL and WCL derivations start from the expan- 
sion of the propagator A for the reduced, system-only 
dynamics, 



p s (t) = A(t,0) Ps (0), 



(•5) 



computed in the interaction picture with respect to 
the renormalized, physical, time-dependent Hamiltonian 
H s (t) =H S + H c {t), where 



H, 



H° s 



(6) 



The renormalizing terms containing powers of A are 
"Lamb-shift" corrections due to the interaction with the 
bath (see, e.g., [38jL The lowest order (Born) approx- 
imation with respect to the coupling constant A yields 
H%° TT , while the higher order terms (• • • ) require going 
beyond the Born approximation. Introducing a cumu- 
lant expansion for the propagator, 



A(t,0)=exp^[A^W(t)], 



(7) 



n=l 



one finds that = 0. The Born approximation con- 
sists of terminating the cumulant expansion at n = 2, 
whence we denote = K: 



A(i,0) = exp[X 2 K (t) + 0(A 3 



(8) 



One obtains 
K(t) Ps 



t <-t 

ds I duF(s — u)S(s)psS(uy 



n 



+ (similar terms) (9) 



as the first term in a cumulant expansion |37| . Here 
F(s) = Tr(pnR(s)R) is the autocorrelation function, 
where pn is the reservoir state and R(s) is R in the Hr- 
interaction picture, and S(u) is S in the interaction pic- 
ture with respect to the physical Hamiltonian Hs(t). The 
"similar terms" in Eq. © are of the form psS(s)S(uy 
and S'(s)S'(u)^ ps- 

At first sight K(t) ~ t 2 , and this is true for small times 
(Zeno effect H^|). The Markov approximation means 
that we can replace K{t) by an expression that is lin- 
ear in t, i.e. 



K(t) 



C(s)ds 



(10) 



where C(t) is a time-dependent Lindblad generator. That 
the Lindblad generator can be time-dependent even after 
transforming back to the Schrodinger picture is impor- 
tant for our considerations below. 



A. Singular Coupling Limit 

The SCL approach we present in this subsection un- 
derlies the standard derivation of the MME that can be 
found in almost any text concerning the Markov approx- 
imation, though not always under the heading "SCL" 
(e.g., [4(|, p.8, Eq. (1.36)). The rigorous derivation of the 
SCL is briefly discussed (with references) in |3J|, pp. 36- 
38. It is based on a rescaling of the bath and system-bath 
Hamiltonians, which physically makes sense in the high- 
temperature limit only. We will shortly see the emergence 
of this limit. 

In essence, the "naive SCL-Markov approximation" is 
obtained by the ansatz F(s) = aS(s) for the autocorrela- 
tion function, whence 



L(s)ps = aS(s)psS(s) 1i + (similar terms). 



(11) 



As a consequence, return to the Schrodinger picture 
gives a MME with the dissipative part independent of 
the Hamiltonian: 



dps 
dt 

£ps 



-i[H s (t),p s ]+Cps, 
~\ 2 a[S, [S, PS ]] 



(12) 



More precisely, we must consider the multi- 
time bath correlation functions F(ti,...,t n ) := 
Tr[p R R(t 1 )...R(t n )} := (R(h)...R(t n )). ' Here 
R(t) := exp(iHnt)Rex-p(—iHfit) are the bath op- 
erators in the interaction picture, pr — exp(—/3Hii)/Z 
(where = 1/kT, Z = Tr[exp(- H R /kT)]) is the bath 
thermal equilibrium state at temperature T, which is a 
stationary state of the reservoir, i.e., [Hji,pfj] = 0. The 
influence of the environment on the system is entirely 
encoded into the {F(ti, ...,t n )}^ =2 . 2 Heuristically, 
the Markov approximation can be justified under the 
following conditions: 

1. The lowest order correlation function, 

/oo 
G{uj)e- iuJt duj, (13) 
-oo 

can be approximated by a Dirac delta function: 3 
F(t) ~ ( / F(s)ds) 5(t) = G(0)5(t) (14) 



(white-noise approximation). Eq. I|13(l defines the 
spectral density G{ui), which is a key object in the 
theory. 



2 F(ti) is constant by stationarity. We reserve the notation F(t) 
for F(ti,t 2 ) = F(t 2 -ti) below. 

3 Note that stationarity implies that F(t) does not depend on s. 
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2. Higher order correlation functions exhibit a 
Gaussian-type behavior, i.e., can be estimated by 
sums of products of the lowest order ones, and then, 
by condition 114fl . decay sufficiently rapidly. 

Let us now comment on the physical relevance of the 
white-noise approximation. 

First, the condition l|14l) cannot be satisfied in general. 
For example, in the important case of linear coupling to 
a bosonic field (e.g., electromagnetic field, phonons in 
solid state), we have G(0) = 0, which means (by inverse 
Fourier transform) that F(t)dt = 0, and therefore 
F{t) cannot be well approximated by 8(t). 

Second, even for models with G(0) > there exists a 
universal relation, the so-called Kubo-Martin-Schwinger 
(KMS) condition, (R(t)R{0)) = (R(0)R(t + i/3)), which 
is valid for all quantum systems at thermal equilibrium. 
This implies: 

G(-uj) = e~ Pu G(uj) . (15) 

(See, e.g., |H[pp.90-91], 41] [pp. 176-177], or H|[p.l37].) 
The fundamental importance of the KMS condition is 
captured by the fact that it is necessary in order for 
thermodynamics to hold. The KMS condition implies 
a strong asymmetry of the spectral density G(u) for low 
T, where T is measured relative to the presence of kT en- 
ergy scales in the bath, i.e., relative to the range where 
G{u>) is non-vanishing. The KMS condition is relevant to 
our discussion since we make the reasonably minimalistic 
assumption that the reservoir (not the QC) is in thermal 
equilibrium. 4 

Third, G(uj) need not be flat even at high T (indeed, 
the KMS condition only implies that G(uS) is symmetric 
at high T). For example, this is the case for the elec- 
tromagnetic field and for phonons, for which at T > 



4 One may challenge the notion that the bath must always be 
in thermal equilibrium. E.g., consider an atom in a microwave 
cavity, with the cavity electromagnetic field initially in thermal 
equilibrium. Now suppose the atom is driven and is coupled to 
the cavity electromagnetic field, which therefore is no longer in 
equilibrium. However, is the internal electromagnetic field the 
relevant environment, or is it the external one? Clearly, the elec- 
tromagnetic field inside the cavity is not a reservoir but itself a 
part of the system. This is because: a) its spectrum is discrete, 
b) its coupling to the atom (close to resonance) is enhanced. 
The reason these considerations matter is because b) implies the 
strong coupling regime, hence failure of the initial state tensor 
product structure assumption, hence difficulties with the sepa- 
ration of the system from the reservoir (dressed atom picture); 
a) implies that F(t) is (quasi)-periodic with short Poincare re- 
currences, hence a strongly non-Markovian regime, and thus as- 
sociated difficulties for Markovian FT-QEC. On the other hand 
the external electromagnetic field has a continuous spectrum and 
the state product structure is easily satisfied, hence qualifies as 
a reservoir. This example merely serves to illustrate accepted 
notions regarding the division into well defined system and bath; 
for most practical purposes a thermal equlibrium is the simplest 
and most relevant model of an environment, and FT-QEC theory 
must be applicable to this setting. 



one has G(u) ~ w 3 /(l - e -^/ feT ) for \u\ < w cut , and 
G(u) — for |cj| > uj cut . One can see that for high T 
(kT 3> ficjeut), G(u>) ~ kTuj 2 is symmetric. Here w cu t is 
the Debye frequency in the case of phonons, while for the 
electromagnetic field w cu t should tend to infinity in the 
renormalization procedure. A flat G(oj) means a struc- 
tureless bath, while physical systems always have a non- 
trivial structure depending on relevant energy scales. 5 

Now let us return to the implications of the SCL as- 
sumptions for the problem of FT-QEC. In order to derive 
the SCL from first principles, one rescales Hr — > iJfj/e 2 , 
rescales Hsr — > Hsu/e, but keeps H$ and pr fixed. 6 
The idea of this rescaling is that it accelerates the reser- 
voir's evolution (via Hr — > Hr/c 2 ) and hence produces 
faster decay of the reservoir correlations, Fit). To see 
this, note that the rescaling Hsr — > Hsr/c increases the 
amplitude F(0) to F(0)/e 2 (proportional to Hg R ), while 
keeping the strength of the noise F(t)dt = G(0) 
fixed (as can be seen via a change of variables t — ► t/e 2 
in the integral). This implies a faster decay of F(t). The 
rescaling procedure is specifically designed to yield the 
delta correlation [Eq. (|14f> ] in the limit as e — > 0. Note 
that if pr is at thermal equilibrium at temperature T 
with respect to Hr, then, since pr — exp(— (3Hr)/Z is 
fixed, it must be at thermal equilibrium with respect to 
Hr/c 2 at the temperature T/e 2 — ► oo, whence our men- 
tion of the high temperature limit, above. Further note 
that Hs is not rescaled since the SCL is (artificially) de- 
signed to produce "white noise" on the natural time scale 
of system's evolution, which is given by Hs- 

Another, equivalent way to understand the emergence 
of the high-T limit is the following: For the Markovian 
condition F(t) ~ aS(t) to hold the spectral density must 
be flat: G(u>) — const. However, this is possible only 
in the limit T — ► oo of the KMS condition. More pre- 
cisely, the Markovian condition can hold only if kT u) 
over the entire spectrum of the system's Bohr frequen- 
cies. Strictly speaking, G(u>) is never constant. The 
variation of G(u>) happens over the "thermal memory" 
time Tth := 1/kT. In the infinite T limit we then re- 
cover the case of zero memory-time, i.e., Markovian dy- 
namics. Physically, it is enough to assume that G(u>) 
is essentially constant over the interval [— ojo,^o] where 
kT > luq ^S> system's Bohr frequencies. I.e., system en- 
ergy scales must be compared to l/Vth and this leads to 
the important realization that the Markovian approxima- 
tion can be consistent with the KMS condition only in the 
high temperature regime kT ^> E, where E is the system 
energy scale. As we argue below, this fact presents a seri- 



5 It is interesting to note that even if we try to enforce a fiat G(uj) 
by, e.g., choosing an appropriate form factor for the spin-boson 
system, the obtained model - the so-called "Ohmic case" - is 
mathematically and physically ill defined (see |42|1. 

6 Note that because different Hamiltonians are rescaled differently, 
this rescaling procedure is not equivalent to a direct rescaling of 
the time variable (which is what is done in the WCL, below). 
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ous difficulty in the context of FT-QEC, the issue being 
essentially that the requirement of a constant supply of 
nearly pure and cold ancillas contradicts the high-T limit 
needed for the Markov approximation to hold. 



B. Weak Coupling Limit 

In the SCL approach above there was no restriction on 
the time-dependence of the system Hamiltonian. How- 
ever, the price paid is the high-T limit. Moreover, while 
mathematically the SCL is rigorous in the scaling limit, 
it is inconsistent with thermodynamics except in the 
T — > oo limit. On the other hand, the derivation by 
Davies, in his seminal 1974 paper [43|. is perhaps the only 
derivation of the MME that is entirely consistent from 
both the mathematical and physical points of view. The 
Davies approach is based not on the high-T limit, but 
rather on the physically plausible idea of weak coupling. 
This is natural and consistent with thermodynamics at 
all temperatures. 

More specifically, Davies' derivation does not invoke 
a flatness condition on G(u) but is, of course, still sub- 
ject to the KMS condition. In the Davies approach the 
Markov approximation is a consequence of weak coupling 
(and hence slow dynamics of the system in the interac- 
tion picture), and time coarse-graining, which leads to 
cancellation of the non-Markovian oscillating terms. The 
price we pay is the invalidity of this approach for time- 
dependent Hamiltonians, except in the adiabatic case. 
We explain this important comment below. Hence, while 
the Davies approach does not require the high-T limit, 
it imposes severe restrictions on the speed of quantum 
gates. 

In his rigorous derivation Davies replaced the heuristic 
condition (|14(l by the weaker 



\F{t)\dt < oo. 



(16) 



This condition avoids the difficulties originating from the 
singularity of the SCL condition (|14f) . and preempts the 
corresponding problems with the high-T limit. 7 ' 8 We now 
consider the cases of a constant, periodic, and arbitrarily 
time-dependent control Hamiltonian. The constant case 



7 In some sense the weak coupling limit is similar to the Central 
Limit Theorem (CLT) in probability, and condition 1161 is anal- 
ogous to a rough upper bound on the second moment in the 
CLT. If it is not satisfied then the noise may be not Gaussian in 
the weak coupling limit. The value of J \F(t)\dt itself does not 
provide any meaningful physical parameter and can depend on 
some regularization/cut-off parameters. 

8 One can go further and ask how generic the Markovian case is, in 
the sense that Eq. 1161 is satisfied. In fact, typically F(t) decays 
as l/t°" (e.g., for the vacuum bath a = 4 .33]), which means that 
in some cases (a < 1) Eq. CJ3 

can be violated. For a systematic 
treatment of these non-Markovian effects see, e.g., 13(1. 



is the one originally treated by Davies 431 , and extended 
in [44! ] to time-dependent Hamiltonians assuming a slow 
("adiabatic") change on the dissipation time scale X 2 t. 
The non-constant cases we study here have, as far as we 
know, not been published before in the general scientific 
literature. 



1. WCL for Constant He: Summary of the Original 
Davies Derivation 

We present a simplified version of the discussion of the 
Markov approximation in j^- Denote by E k the Bohr 
energies (eigenvalues of H s ), let u E {uj k i = E k - Ei} k ,i, 
and let S u be the discrete Fourier components of the 
interaction picture S, i.e., 

S(t) = exp(iHst)S exp(— iiJsi) = 5 W exp(icji), 



(17) 

where Hs is the renormalized (physical) system Hamil- 
tonian: the sum of the "bare" H s and a Lamb shift term 
(bath induced), as in Eq. ©. Equivalcntly, 

[Hs,S u ]=vS u . (18) 

We remark that in the original Davies paper the Bohr 
energies and Eq. (|18|l are computed with respect to the 
bare Hamiltonian Hg. Here we use the physical Hamil- 
tonian Hs in order to take into account the fact that the 
Lamb shift term, although formally proportional to A 2 , 
can be large or even infinite after cut-off removal. 
Then, it follows from Eq. (JSJ that 

K(t) Ps = y^SupsSl [ e^-^du [ " F{T)e^ T dT 

+ (similar terms). (19) 

The weak coupling limit is next formally introduced by 
rescaling the time t to t/X 2 (van Hove limit). This en- 
ables two crucial approximations, which are valid in the 
resulting large- t limit: 

1. We replace 9 

This makes sense for 

t > max{l/(w - a/)}. 



(20) 



(21) 



This violates Al, expressed in terms of the Bohr 
frequencies. We see here already the emergence of 
an adiabatic criterion for the validity of the Markov 
approximation. 



be used, but the result is essentially the same 



ncipa 
37]. 
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2. We replace J_ " F(r)e lWT dr by the Fourier 
transform: 

/t—u poo 
F{r)e lu)T dT r> = / F{r)e lu)T dT. (22) 

The physical validity of the last approximation is usu- 
ally ignored, though one can make the following argu- 
ment: On the LHS of Eq. (|22fl . for a given Bohr frequency 
u the Fourier-like integral must sample the function F(t) 
with sufficiently high accuracy so that the Fourier trans- 
form approximation will be valid. To this end one needs 
a time t such that: (i) t 3> l/u>. This is a weaker 
condition than the previous one [t 3> max{l/(o; — w')}] 
which involves differences of Bohr frequencies, (ii) The 
time t must be also much longer than the time scale 
of the wildest variations of F(t), which is typically [as 
may be checked for simple models of spectral densities 
G(u>)] given by l/w cut , where aj cut is a high-frequency 
cutoff. When lu < uj cut (i) implies (ii). Therefore typi- 
cally Eq. I|2(J|) is a stronger assumption than Eq. I|22|) . 

Applying the approximations (|20() and (|22|l . we obtain 
K{t)ps = t^2 u S^PsSIjG(uj) + (similar terms), and hence 
it follows from Eq. (|f 0|l that C(s) = C is the Davies 
generator in the familiar Lindblad form: 

— T- = ~i[Hs, ps] + £ps, 
dt 

C PS = Ix^G^dSu^SD + iSuP^l}) (23) 
Several remarks are in order: 

(i) The absence of off-diagonal terms in Eq. com- 
pared to Eq. I jl 9|) . is the hallmark of the Marko- 
vian limit. Namely, the Davies derivation relies 
on the cancellation of the non-Markovian off-diagonal 
terms £ wA ,, SupsS^ f e^-^du. This time coarse- 
graining is possible due to integration over the fast- 
oscillating f Q e l ^ u ~ u terms over a long timescale, i.e., 
over t ^S> max{l/(w — ui')} (see also 38]). As remarked 
above, this violates Al, expressed in terms of the Bohr 
frequencies. 

(ii) It follows from Bochner's theorem applied to the 
Fourier transform definition of G(uj) that G{u>) > 
H3[P-90], (3||[p-136]; this result is essential for the com- 
plete positivity of the Markovian master equation in the 
WCL. 

(iii) Davies' derivation showed implicitly that the notion 
of "bath's correlation time" is not well-defined - Marko- 
vian behavior involves a rather complicated cooperation 
between system and bath dynamics. More specifically, 
the relations and 118|) together imply that the noise 
and Hs are strongly correlated. In other words, contrary 
to what is often done in phenomenological treatments, 
one cannot combine arbitrary Hs 's with given S u 's. This 
point is particularly relevant in the context of FT-QEC, 
where it is common to assume Markovian dynamics and 
apply arbitrary control Hamiltonians. 



Davies did not consider time-dependent system Hamil- 
tonians in [43j], but it is possible to generalize his deriva- 
tion to allow for slowly varying system Hamiltonians 
[37L [51 145| . That is, whenever the time scale of the vari- 
ation of He {t) is much longer than the inverse of the typ- 
ical Bohr frequency (of Hs), it is possible to add Hc(t) 
to the system Hamiltonian in Eq. I|23|) . necessitating at 
the same time this change also in Eq. I)18|l . This is a type 
of adiabatic limit (indeed, the in Eq. (fT5|) can be in- 
terpreted, with Hs replaced by Hs + Hc(t), as being adi- 
abatic eigenvectors of the superoperator [Hs + Hc(t), •]). 
We note that an alternative approach to adiabaticity 
in open quantum systems was recently developed in 
Ref. ^(J- This approach, while being very general, is 
more phenomenological in that it postulates a convolu- 
tionless master equation, and then derives corresponding 
adiabaticity conditions. Closer in spirit to the Davies 
derivation is another recent approach to adiabaticity in 
open systems, which assumes slow system variation to- 
gether with weak system- bath coupling pTf. 



2. WCL for Periodic Driving: Floquet Analysis 

Before considering the case of periodic He let us con- 
sider briefly once more the case of a constant Hamilto- 
nian in the so-called covariant dissipation setting. Co- 
variance is the commutation condition TLC = CTL where 
H = [Hs, ■} is the super-operator constant Hamiltonian, 
and C is the Davies generator. Covariance is an abstract 
property which is automatically fulfilled for the Davies 
generator. 10 It is convenient since it implies factoriza- 
tion of the full propagator into Hamiltonian and dissipa- 
tive parts. Markovian dynamics obtained in the WCL as 
discussed above takes the form 



-^ = {-iH + C)p, t>0, (24) 

where the most general form of the Lindblad (or Davies) 
C satisfying Eq. (J2U is 

Here {ui} = Spectrum(7i), i.e., the Bohr frequencies (dif- 
ferences of eigenvalues of H), and 

W j («)=wV j (w) (26) 



This can be verified by directly computing "HC and making use 
of Eq. (SJ and the relation [A, BC] = [A,B]C + B[A,C] (for 
operators A, B and C). A more elegant way to see this is to 
consider C(t) = exp(— ifH)C exp(ifH) and note that Eq. 1181 
implies that S(t) and (t) rotate in opposite directions. Hence 
C(t) = C, whence dC(t)/dt = gives the result. 
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[i.e., Eq. 1|18|) ]. The solution, i.e., the dynamical semi- 
group is 

p(t)=A(t)p(0), A(t) = e- ltn e tc . (27) 
Now consider a periodic control Hamiltonian with period 

e 

H c (t) =Hc(t + 0), Q = 2n/e. (28) 

(Note that f2 is not the Rabi frequency, which throughout 
this paper we denote by SIr.) The situation is then very 
similar to the standard (time-independent He) WCL, 
but the set of "effective Bohr frequencies" (Floquet spec- 
trum) lu is now larger and is of the form {lu + q£l}, 
q = 0, ±1, .... Here lu are Bohr frequencies for the Flo- 
quet unitary [defined in Eq. ll.'U)[l below], i.e., differences 
of eigenvalues e a of the Floquet unitary, rather than 
{lu} = Spectrum(7i) as above. As this set of "effec- 
tive Bohr frequencies" is discrete the WCL still works, 
but the final Davies generator is more complicated, as 
we now show. 

Define the time-ordered unitary propagator 

U(t, s) = T exp (-i jf H s (u)du^j , t>s (29) 

which satisfies the properties U(s,t) = U(t, s) _1 = 

u{t, s )\ u(t, s )u(s,u) = u(t, u ), u(t + e,s + e) = 

U(t,s), and f t U(t,s) = -iH s (t)U(t, s), f t U(t,s)^ = 
iU(t, s)^ Hs(t). The Floquet unitary operator is 

F(s) = U(s + Q,s), (30) 

with corresponding super-operator action 

T(8)p = F{a)pF(a)\ (31) 

and Floquet eigenvectors \<f> a ) and eigenvalues (quasi- 
energies) e a satisfying 11 

F(Q)\<f> a ) =e-^ e|0 Q ). (32) 

It follows from standard Floquet theory that 

U(t,0)\</> a ) = z~ ltta E e~ Uqn \M<l)), (33) 
qez 

i.e., the set {\(/>a(q)}} is a complete basis. Therefore we 
have at most as many g's as the dimension of the Hilbert 
space. That the number of g's is finite is important for 
our considerations below. 



11 Note that the Floquet Hamiltonian H,s(t) — id/dt operates on a 
different Hilbert space than F(Q) (the space of periodic functions 
with values in the system's Hilbert space). But its eigenvalues 
coincide with t a from Eq. 1321 . 



We call a Lindblad generator L a "covariant dissipative 
perturbation of Hs(t)" if 

T{Q)C = £JF(0) (34) 

We will assume this property, similarly to the case of 
a constant Hamiltonian described above. In fact, covari- 
ance holds for a periodic Hs(t) and also for the corre- 
sponding WCL Davies generator. One can then derive 
the covariant master equation (we sketch the derivation 
below) : 

J = {-iH{t) + £(*)) p, t>0, (35) 

[compare to Eq. (|24|l ] where 

C(t) = U{t,0)CU(t,0) f , (36) 
■jM{t,s) = -iH(t)U(t,s), (37) 

and where the general form of C appearing in Eq. (|36|l is 
given by Eq. I|25|l . with Vj(lu) now being eigenvectors of 

F(O)V j (cj) = e- i " V j (Lj), (38) 

rather than of Tt, as in Eq. I|26(l . Moreover, here 
{lu} = {e a — where e Q are quasi-energies (effec- 
tive Bohr frequencies) of the Floquet operator, rather 
{lu} = Spectrum(7i) as we saw in the case of constant 

The solution replacing Eq. (|27|l is 

pit) = A(t,s)p(s), t>s 
A(t,s) = Texpj^ {-iH{u) + £(u)) duX (39) 

By direct computation one can prove the following prop- 
erties: 



C(t + Q) 


= C(t), 


(40) 


^(s)C(s)^{s) ii 


= C(s), 


(41) 


A(t, s)A(s,u) 


= A(t, u) for t > s > u, 


(42) 


A(t + Q,s + Q) 


= A(t,s), 


(43) 


A(t,s) 


= U{t,s)e- {t - s)c{s] . 


(44) 



To derive the covariant master equation (|35|l one con- 
siders the standard picture of an open system S + R with 
the total Hamiltonian 

H SR (t) =H° s {t)+H R + Y,Sk®Rk, (45) 

fe 

(we neglect the Lamb shift correction here; it can be 
included, changing Hg(t) into the physical Hamiltonian 
Hg(t), by a suitable renormalization procedure), station- 
ary reservoir state p R , [H R ,p R ] = 0, Tr(p R -) = (•)#, 
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(Rk)n = 0. Then, exactly following a Davies-like calcu- 
lation using a Fourier decomposition of S(t), now gov- 
erned by a periodic Hamiltonian, and making in partic- 
ular again the crucial assumption Eq. 121|) . which now 
reads 

t»max{l/(ct>-a/-r-mO)}, m = 0, ±1, ±2, ... (46) 

with \m\ upper-bounded by the dimension of the Hilbcrt 
space [see remark after Eq. (|33|l ]. one obtains Eq. ilool) in 
the Davies WCL. The explicit form of the generator is: 

k,l qeZ {w} 

+ w),pSfc(<?, (47) 

Here {cu} = {e a ~ ep}, the Floquet spectrum, and 

/oo 
e- ttx (R k (t)Ri) R dt (48) 
-OO 

S k (q,ui) = 'J2 ^2 (<p a (p + q)\S k \(f> a >{p))\(f> a )((t> a >\- 

p6Z {e a -e a ,=uj} 

(49) 

Sk{q,io) is the part of S(t) which rotates with frequency 
cu + q^l and can be computed using Eq. J33J. Note that by 
diagonalizing the matrices Rki one can transform the gen- 
erator C of Eq. i|T7|) into the form of Eq. , which al- 
lows one to read off the operators Vj(co) appearing there. 
Now to some important comments: 

- Timescale analysis: Note that for the periodic case the 
differences of "Bohr frequencies" may be of the order of 
1/0. Hence we conclude from Eq. (|46[) that one must 
average over many periods O, i.e., require f > 8. This 
can be interpreted as a condition that "the environment 
must learn that the Hamiltonian is periodic" . This is 
exactly analogous to the adiabaticity condition in the 
adiabatic case: H(t) must be constant over many inverse 
Bohr frequencies to "be recognised" by the environ- 
ment. The periodic WCL is also a coarse-grained time 
description with the additional time scale 0. Note that 
arbitrarily fast periodic driving ( small Q) is incompatible 
even with the kind of generalized, finitely localized MME 
derived here, since then differences of Bohr frequencies 
matter in Eq. \4°\ l (recall that max|m| is bounded by 
the - typically small - dimension of the system Hilbert 
space) . 

- Where is the Rabi frequency? Note the dependence of 
the operators Sk(q,to) on the Floquet eigenvalue differ- 
ences e a — e a i . The usual Rabi frequency, fin = 2dE/h 
(d is the dipole moment, E is the electric field ampli- 
tude) arises in the dipole approximation, which we have 
not made here. The usual Rabi frequency is replaced in 
our non-perturbative treatment (in the sense of no multi- 
pole expansion) by the difference of Floquet eigenvalues 



e a - e a , in Eq. JUJ. 12 

- More on the Rabi frequency: As we saw, the non- 
Markovian terms vanish because of the time coarse- 
grained description. To attain this, we must average over 
times t 3> max C j ^/{l/(aj — ^>')}, but must also keep in 
mind that the longest relevant scale for coarse-graining is 
given by the exponential decay time r (a derived quan- 
tity), i.e., we must have t < r. The Rabi frequency Or is 
a difference of two Bohr frequencies u,u': Or = ui — uj' . 
This implies that coarse-graining does not makes sense 
if Ort <§; 1 [since then t < t <C I/Or = l/(w — u'), in 
contradiction to the fundamental requirement on t]. In 
physical terms this means that the width of the spectral 
line (7 = 1/t) is larger than the level splitting Or (see, 
e.g., Fig. 2.5 (i),(h) in |4jj f° r an illustration in the case 
of the incoherent fluorescence spectrum) and therefore 
"the environment has no time to recognize the details 
of the spectrum" . On the other hand, when Ort 3> 1 
(not inconsistent with the WCL) , Or must appear in the 
generator, as appears from our treatment of the case of 
periodic driving in the WCL, above. Unfortunately there 
are examples in the literature where an MME is written 
down subject to Ort ^> 1 but Or does not appear in the 
generator [e.g., Eq. (2.96) in where Or ~ 10 10 Hz 
and t ~ lCT 8 s]. 

- Quantum optics considerations: The Markov approxi- 
mation is commonly accepted as an excellent approxima- 
tion in quantum optics; see, e.g., the discussion of res- 
onance fluorescence in [40j[Ch.2]. This is also the basis 
for substantial confidence in the possibility of FT-QEC in 
quantum optical systems, such as trapped ions [48j and 
atoms trapped in microwave cavities |49j. Such argu- 
ments are based on the relative flatness of the damping 
constants j(u>) as a function of frequency. This argu- 
ment is closely related to the notion of the flatness of 
the spectral density G(lo) in the SCL, since the damping 
constants are proportional to G(ui) [see Eq. For 
example, below Eq. (2.95) in ^(j the author argues that 
one can write down a Rabi frequency-independent MME 
for resonance fluorescence since 7(1^,4) and j(uja ± Or) 
(where to a is the Bohr frequency) differ by less than 
0.01% at optical frequencies and reasonable laser inten- 
sities. However, this ignores the corrections due to the 
Rabi frequency to the operators Sk.(q,uj) [Eq. l|4*9"|) ]. This 
disagreement can be traced to the question of at which 
point in the derivation it is safe to neglect Or; in 
this is done on the basis of the flatness of 7(0*) before "a 
lot of tedious algebra" f4(j[p.48], but our Floquet analy- 
sis shows that, in fact, one cannot neglect the Rabi fre- 



One can see that such a term also arises in th e us ual dipole ap- 
proximation by considering, e.g., Eq. (2.94) in |4(1 . The interac- 
tion picture raising and lowering operators cr± (t) (for a two-level 
atom driven by a classical field) there oscillate with three "Bohr 
frequencies" u>a,u>a ± , where Q = 2dE/h denotes the usual 
Rabi frequency. Hence the Rabi frequency is a difference of two 
Bohr frequencies. 
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quency relative to the Bohr frequency. This is relevant for 
our general discussion since the "finitely localized MME" 
which is the outcome of the Floquet analysis (see next 
comment) actually exhibits a weak non-Markovian char- 
acter. Such deviations are, of course, important for FT- 
QEC, even if the effects are small. We revisit this point 
in Section Ivl below. 

- Are there any non-Markovian effects at work here? It 
seems that one should accept the generalized notion of a 
quantum Markovian master equation as the one given by 
Eqs. lj55T) . and J2HJ), i.e., a master equation with a 
possibly time-dependent Lindblad generator. In Davies' 
generalization to the time-dependent case 0] ("adia- 
batic WCL") the dissipative generator C{t) depends on 
the Hamiltonian at the same time t. This is a type of 
"local generalized MME". On the other hand, in the pe- 
riodic WCL treated here, the dissipative generator C(t) 
depends on the Hamiltonians Hs{u) from an interval, 
say [0,t] (t < 8), as can be seen from Eq. (|36[) . which 
involves U(t,0). This is therefore a type of "finitely lo- 
calized MME" , though one could argue that it exhibits a 
weakly non-Markovian character because of this depen- 
dence of the dissipative generator on the past. On the 
other hand, a non-Markovian master equation (in the 
convolutionless formalism j^) is also given by Eq. 
but the generator is not of Lindblad form [in particu- 
lar, it is not of the form l|32f> ]. and may depend on the 
Hamiltonian in the distant past. The weight of distant 
past contributions depends on the decay properties of 
F(t) which are, generically, not exponential but rather 
powerlike. In the WCL the non-Lindbladian terms van- 
ish due to the oscillating character of the e 4 '" - " ^ terms 
in Eq. $P$ . 

- The original Davies derivation: We note that the 
Davies result is a limit theorem which states that for a 
sufficiently small coupling constant the WCL semigroup 
is a good approximation to the real dynamics. However, 
Davies' theorem itself does not provide the conditions 
under which a given physical coupling is "small enough" . 
In particular, one cannot extract from Davies' theorem 
under what conditions the fast oscillating terms vanish. 
This can, however, be done by a more heuristic analysis, 
as done above. 



3. WCL for an Arbitrary Pulse 

We now consider the case 

H c (t) = H + f(t)H 1 , (50) 

i.e., an arbitrary driving field. This is, of course, the case 
of most interest in FT-QEC. It follows from Fourier anal- 
ysis that this case can be treated qualitatively as a "su- 
perposition" of periodic perturbations discussed above. 
For a single frequency f2 the validity of the Marko- 
vian approximation is restricted by the condition 1)46(1: 
t max{l/(oj — lo' + mfl)}. The discreteness of the 
frequencies {lo} and {mfi} is key: it allows for condition 



(|46ll to be satisfied with finite t. A pulse f(t) has a con- 
tinuous band of frequencies of width T ~ l/r g (where 
T g is the gate duration), with amplitudes (Fourier trans- 
form) /(fi), which add to and smear the effective Bohr 
spectrum {lo}. If the pulse is long (a slow gate) then only 
a narrow band appears, and the smearing effect is unim- 
portant. More precisely, if l/r g is much smaller than 
the typical difference of the Bohr frequencies, the "en- 
ergy quanta" mfi [with |m| restricted by the (typically 
small) dimension of the system Hilbert space] cannot fill 
the gap between lo and lo' and the condition 14tj|) can be 
satisfied. This is our adiabatic approximation. For fast 
pulses, when l/r g is comparable to \lo — lo'\, the condi- 
tion H4tjfl cannot be fulfilled: the effective Bohr spectrum 
becomes quasi-continuous and the denominator in condi- 
tion 14(j|) becomes abitrarily small. The result is that the 
WCL analysis breaks down and non-Markovian effects 
dominate. 

Thus, the condition for the adiabatic limit (Markov 
approximation valid) is: "the width of the band is much 
smaller than the minimal difference of the effective Bohr 
frequencies" . This is in contradiction with the fast gate 
assumption, Al. 



C. Section Summary 

The main advantage of the MME (|23|l is its consis- 
tency with thermodynamics. Namely, as a consequence 
of the KMS condition l|15|) and the condition l|18|) . for 
a generic initial state the system tends to its thermal 
equilibrium (Gibbs) state at the temperature of the heat 
bath [35| . (An important exception to this rule are states 
within a decoherence-free subspace E], but these 
states are not generic due to required symmetry proper- 
ties of the system-bath interaction.) Therefore the dis- 
sipative part of the generator must depend strongly on 
the Hamiltonian dynamics. This is consistent with the 
notion of a coarse-grained description familiar from the 
study of MMEs: the bath needs a time much longer than 
max Utl 1/lom to "learn" the system's Hamiltonian in or- 
der to drive it to a proper Gibbs state. In other words, the 
Markov approximation is, equivalently, a long-time limit 
(compared to max Ull 1/lom - the system's Bohr frequen- 
cies), and one cannot expect this approximation to be 
valid at short times. However, FT-QEC assumes oper- 
ations on a time-scale that is short on the scale set by 
max Wfei 1/oJm- 

Strictly speaking the MME is valid only when Hs 
is not time dependent. As we have shown, one can re- 
lax this by assuming slowly varying Hs, giving rise to 
an "adiabatic MME", Eqs. and How- 

ever, to accept Eqs. (|37)|l and as a genuine 

Markovian description is somewhat of a stretch, since 
the real question is not whether one obtains the Lindblad 
form, but rather how C(t) depends on the Hamiltonians 
Hs(u), locally (i.e. u ~ t) or nonlocally. For fast gates 
and generic environments the dependence is non-local, 
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involving memory effects. In any case, the crucial con- 
dition that must be satisfied for a (generalized) MME is 
Eq. H46(l . which implies that the average Bohr spectrum 
must be discrete. In essence, as long as the applied con- 
trol does not spoil this discreteness a (generalized) MME 
can be derived. On the other hand, this means that fast 
gates are incompatible with the MME, in violation of Al 
of FT-QEC theory. The corollary: finite speed of gates 
implies non-Markovian effects. 



IV. ARE THE STANDARD FT-QEC 
ASSUMPTIONS INTERNALLY CONSISTENT? 



We now briefly summarize our examination of the as- 
sumptions of FT-QEC, in light of the considerations 
above, and highlight where there may be internal incon- 
sistencies in FT-QEC. As discussed above, there are es- 
sentially two rigorous approaches to the derivation of the 
MME: (i) the SCL, which is compatible with arbitrarily 
fast Hamiltonian manipulations, but requires the high-T 
limit; (ii) the WCL, which is compatible with thermody- 
namics at arbitrary T, but requires adiabatic Hamilto- 
nian manipulations. 

The standard theory of FT-QEC (excluding Refs. [H 
H3>E1) requires a quantum computer (QC) undergoing 
Markovian dynamics, supplemented with a constant sup- 
ply of cold and fresh ancillas. These assumptions are 
contradictory under the SCL, since the QC would have 
to be at high-T, while the ancillas require low-T on the 
same energy scale E (set by the Bohr energies of the sys- 
tem = computer + ancillas). Specifically, if we were to 
assume that for the ancillas too kT ^S> E, they would 
quickly become highly mixed. If we insist that kT <C E 
for the ancillas, then by coupling them to the QC we can 
no longer assume, in the SCL, that the total system = 
QC + ancillas is described by Markovian dynamics. 

If, on the other hand, we approach the problem from 
the (physically more consistent) WCL, then A3 is in- 
compatible with Al (the assumption of fast gates). 
Namely, in the WCL only adiabatic Hamiltonian manip- 
ulations are allowed. Specifically, the Markov approxi- 
mation in the WCL requires a discrete system (effective) 
Bohr frequency spectrum, such that the condition r g 3> 
max,j tl 1/wfei can be satisfied, hence violating the t 9 lob = 
0(tt) condition of Al. These conclusions are unavoid- 
able if one accepts thermodynamics, since they follow 
from seeking a Markovian master equation that satisfies 
the KMS condition - a necessary condition for return 
to thermodynamic equilibrium in the absence of external 
driving. We take here the reasonable position that a fault 
tolerant QC cannot be in violation of thermodynamics. 



V. POSSIBLE OBJECTIONS TO THE 
INCONSISTENCY 

In this section we analyze a list of possible objections 
to the inconsistency we have pointed out. 



A. Is thermodynamics relevant? 

With respect to the SCL: "Thermodynamics is irrel- 
evant (since a QC need not ever be in thermal equilib- 
rium)." 

Note that we never claim that the QC is in thermal 
equilibrium; only the bath is. This assumption is a sim- 
plification which allows us to use a single parameter T 
and therefore a single "thermal memory time" h/kT. 
There is no reason to use a nonthermal bath or many heat 
baths with different temperatures: this does not make the 
spectral density flat and can only introduce more param- 
eters. 



B. Doesn't the interaction picture save the day? 

With respect to the WCL: "Suppose we have the fol- 
lowing Hamiltonian in the Schrodinger picture: H = 
H s + H c {t) + H S R + H R where \\H S \\ > \\H C \\ =control 
Hamiltonian \\Hsr\\- Then in the interaction pic- 
ture with respect to H$ the term He is dominant and 
hence can implement fast gates. However, in the original 
Schrodinger picture He is small and hence the adiabatic 
limit for the derivation of the MME is satisfied. Thus 
we have an example where we can have fast gates ( in the 
interaction picture) and still the WCL can be satisfied so 
that the Markovian limit can be reached. Moreover, this 
is the relevant limit relevant for quantum optics, e.g., 
trapped ions." 

There are a number of problems with this argument. 
First, one should be more careful about the formulation 
of the condition for adiabaticity. It can be stated as 
\du(t)/dt\ <C oj(t) 2 , where u)(t) is a "relevant" Bohr fre- 
quency. Merely comparing norms as above does not guar- 
antee adiabaticity. Second, in the quantum optics con- 
text we note the following. For three-level trapped ions 
we have two Bohr frequencies: a large, time-independent 
u>i, and a small, time-dependent W2(t) (degenerate lev- 
els splitting). Only 0J2 is "relevant" because it is re- 
lated to gates, and then the adiabatic condition implies 
that \du)2{t)/dt\ is correspondingly small, which contra- 
dicts the fast gate condition Al. Third, the inequality 
1 1 flic 1 1 ||flsfl|| is in fact not satisfied in the Markovian 
WCL, where ||flsfl|| diverges (one should not confuse 
the small system-reservoir coupling parameter involved 
in the van- Hove limit with the operator norm, which can 
be infinite). 
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C. Doesn't quantum optics provide a 
counterexample? 

With respect to the WCL: 11 Trapped ions and other 
quantum optics systems provide a counter-example: a 
system experimentally satisfying Markovian dynamics 
and allowing fast Rabi operations." 

We have already addressed quantum optical systems 
in Section IIII B 21 Let us add a few comments. We 
do not know of any quantum optics experiment test- 
ing the Markov approximation with the accuracy rele- 
vant for FT-QEC (for quantum dots, on the other hand, 
non-Markovian effects are very visible). We know that 
for constant, and also for strictly periodic Hamiltonians 
(which corresponds in quantum optics to a constant ex- 
ternal laser field), the Davies derivation can be applied 
(or extended, as in Section IIII Bl) and the Markov ap- 
proximation is applicable. The problem appears for fast 
gates. It would be difficult to test the Markov approxi- 
mation in this case with the required accuracy, because, 
e.g., the results depend on the shape of the pulse. A rel- 
evant example is resonance fluorescence, as described in 
H3[pp.43-61], and as discussed in Section HITrT2l The 
damping effects are only present in the widths of spec- 
tral lines - see ^3[p.61, Fig- 2.5]. The Markov approxi- 
mation gives Lorentzians while non-Markovian dynamics 
may give rise to more complicated lineshapes. Consider 
a 2-level atom like in |4fJ Section 2.3.2., and in partic- 
ular the final formula Eq. (2.96), which describes reso- 
nance fluorescence via a MME. The author claims that 
for typical parameters in quantum optics the dissipative 
part does not depend on the Rabi frequency J7r [recall 
our discussion in Section fill B 2) . Hence, as the gates 
are entirely related to Or, it appears that either fast or 
slow gates are possible. The argument is based on the 
small ratio Or/w^ < 10 10 /10 15 (where uja is the Bohr 
frequency). This is fine for replacing the spectral density 
at uja ± by the density at uja, but the subsequent 
argument that we can replace [in Eq. (2.94)] £1r by is 
inaccurate. This would be correct only if the decay time 
t = I/7 is short enough such that Qrt <C 1. However, as 
explained in Section Till B 21 in this case the Davies type 
averaging makes no sense physically. In fact, typically 
for radiation damping r = 10 _8 s, and then SIrt < 100 
only. Hence for a fixed Or we do in fact not have a sim- 
ple Lindblad generator (of the type (2.96) in but 
rather a more complicated generator with Lindblad op- 
erators depending on the Rabi frequency, as in Eq. (|4"7Jl . 
Again, in the derivation of a proper generator an aver- 
aging over terms of the form exp(— iQnt) must be per- 
formed. Therefore the condition for the adiabatic ap- 
proximation involves the Rabi frequency Or and cannot 
be satisfied for fast gates. For experiments based on spec- 
tral measurements the difference between the two types 
of generators we have just discussed is probably irrel- 
evant for many reasons; however, the quantum state of 
the atom at a given moment is sensitive to a small change 
in the Lindblad operators, and this is important in a fault 



tolerant implementation of quantum logic gates. 



D. Is Al truly an assumption of FT-QEC? 

With respect to the WCL: "Doesn't Al impose an un- 
necessary constraint on FT-QEC, in that gates are not 
required to satisfy the condition T g u> = O(ir) ?" 

In other words, one might argue in favor of slow gates, 
where instead the condition is t 9 uj 3> 0(ir). Such gates 
are certainly relevant in the context of the adiabatic 
quantum computing (AQC) paradigm psj, holonomic 
QC [53I l54l . or topological quantum computing (TQC) 
pai l5a . 157| . We comment in more detail on AQC, HQC, 
and TQC in Section Ivl! The question of interest to us is 
whether an adiabatic gate satisfying t 9 lu 3> 0(ir) is appli- 
cable to the standard FT-QEC paradigm we are consid- 
ering here, and which is very different from AQC, HQC, 
and TQC. 

First, let us clarify that by gates we mean one and 
two-qubit unitaries picked from well-known discrete and 
small sets of universal gates j5g. An algorithm is con- 
structed via a sequence of such gates, and computational 
complexity is measured in terms of the minimal num- 
ber of required gates. Of course one can instead join all 
gates used in a given algorithm into a single unitary and 
call this a gate, but then one runs into the problem of 
finding a relevant (physical) Hamiltonian and quantify- 
ing computational complexity. For a given gate there are 
infinitely many Hamiltonian realizations. Among these 
are fast ones (optimal) which satisfy t 9 lu — 0(tt) and 
slow ones (adiabatic) satisfying r g uj ^S> 0(ir) (all inequal- 
ities here are in the sense of orders of magnitude). For 
example, consider a 7r-rotation. The fast (optimal) real- 
ization satisfies r g uj — tt (compatible with Al), while the 
slow (adiabatic) one satisfies r g uj = tt + 2im with n»l 
(contradicts Al). 

Now, one may ask whether a slow realization of gates 
can prevent the inconsistency with the WCL. We argue, 
based on computational complexity considerations, that 
the answer to this question is negative. To see this, note 
first that non-Markovian errors are uncorrectable in stan- 
dard FT-QEC. Therefore such non-Markovian, uncor- 
rectable errors accumulate during the computation (by 
definition, they are not corrected by "Markovian FT- 
QEC"), and in order to keep them under control, the 
probability of such errors per gate, Pnon-M, should scale 
as 

Pnon-M ~ 0[1/ (volume of algorithm)] 

= 0[l/(input size)"], (51) 

where a is some fixed power. Now, it follows from our dis- 
cussion m Section llTTBl that the more adiabatic the evolu- 
tion, the smaller is the probability of the non-Markovian 
errors per gate. Therefore, if one writes the adiabatic- 
ity condition as r g uj > M, where M ^> 1 is the "adia- 
batic slowness parameter" , then the probability of non- 
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Markovian errors should satisfy 

0{l/M?), (52) 

where (3 is another fixed power [lu (the Bohr or Rabi 
frequency) is limited essentially by the choice of physical 
system]. Comparing the two expressions for p n on-M, we 
see that M must grow with input size. This means that 
if one works with adiabatic gates in order to keep the 
dynamics (approximately) Markovian, the result is that 
one must slow the gates in proportion to the input size 
(to some power). This, however, violates the threshold 
condition of FT-QEC, in which the input size and gate 
times are independent parameters (see, e.g., Theorem 12 
in 0). 



E. Measurements 

With respect to both the WCL and the SCL: "Recent 
results on fault-tolerant QC using measurements only 
(e.g., JEUi wdtj ) render all the claimed problems irrele- 
vant." 

Indeed, we have so far discussed only the problems with 
quantum logic gates. Moreover, measurements are an 
integral part of FT-QEC theory as well, in particular to 
reset and disentangle ancillas before they are introduced 
into an error-correction circuit. Therefore some remarks 
on the use of measurements are in order. 

In the most advanced FT-QEC scheme of 0, mea- 
surements are performed at the end of the computation. 
However, this approach demands a high resource over- 
head, which may make it impractical. Therefore, more 
recent proposals (e.g., [IjllsH) Te ^f on feedback mecha- 
nisms employing the results of quantum measurements. 
Those "measurements in the middle of computation" are 
treated for simplicity as certain von-Neumann projective 
measurements (but with efficiency < 1) satisfying a re- 
peatability condition. The latter implies that the subse- 
quent measurements reduce the measurement error ex- 
ponentially as their number increases. This assumption 
should be carefully scrutinized, within realistic Hamilto- 
nian models of quantum measurement treated as a dy- 
namical process. Here, again one can expect that the 
tacit assumption of statistical independence of repeated 
measurements is in conflict with the non-Markovian char- 
acter of the dynamics of open quantum systems. 

As all proposed measurement schemes are based on 
electromagnetic interactions, it should be possible to con- 
struct a rather general Hamiltonian framework and ap- 
ply it to various particular implementations. Indeed, this 
has been done, e.g., for a single-electron tunneling (SET) 
transistor coupled capacitively to a Josephson junction 
qubit 62] . Rather than assuming that the measurement 
apparatus is coupled to the system whenever measure- 
ments must be performed - an option which is hard to 
achieve in mesoscopic systems - Ref. [62j ] makes the rea- 
sonable assumption that the measurement apparatus is 



always coupled to the system, but is in a state of equi- 
librium when it is not needed. A measurement is then 
performed by driving the measuring device out of equi- 
librium, in a manner that dephases the qubit to be mea- 
sured. Generic features emerging from this analysis are 
the existence of three different time-scales characterizing 
the measurement: the dephasing time, the measurement 
time (which may be longer than the dephasing time), 
and the mixing time (the time after which all the in- 
formation about the initial quantum state is lost due to 
the transitions induced by the measurement). Ref. |62l ] 
thus arrives at a criterion for a "good" quantum measure- 
ment: the mixing time should be longer than the mea- 
surement time. A time-scale analysis of measurements 
in optical systems, accounting for spontaneous emission, 
can be found, e.g., in Ref. [63j. A fully consistent anal- 
ysis of FT-QEC should account for the existence of such 
time-scales in a dynamic description of the measurement 
process. In particular, it is important to set appropriate 
bounds on these time-scales, so that they may be taken 
into account in a threshold calculation (an analysis based 
on a stochastic error model was reported in Ref. |!4jL 



F. Degenerate Qubits 

With respect to the SCL: "Degenerate qubits automat- 
ically satisfy the high T limit since their intrinsic energy 
scale vanishes." 

Examples of degenerate qubits are common, e.g., in 
trapped ion quantum computing implementations where 
a pair of degenerate hyperfine states can serve as a qubit, 
with an auxiliary third level used to implement quantum 
logic gates via Raman transitions [64[ • The case of degen- 
erate qubits is somewhat more subtle to analyze within 
the context we have explained above. Naively, in such 
a case the high-T limit is indeed automatically satisfied, 
since the system energy scale is zero. Therefore it appears 
that one could claim that the SCL version of the Markov 
approximation is attainable. However, upon closer ex- 
amination this still seems problematic. Indeed, the van- 
ishing of an energy scale for degenerate qubits holds, 
strictly speaking, only for fully adiabatic techniques, e.g., 
HQC |53t l54|. Otherwise transformations between logi- 
cal states are achieved by resorting to effective Hamil- 
tonians which involve virtual transitions. For instance, 
if |0) and |1) denote degenerate qubit levels (e.g., hy- 
perfine levels of an ion), one can introduce far-detuned 
(e.g., laser) couplings of |0) and |1) with a third auxiliary 
level. Second order perturbation theory then yields the 
effective Hamiltonian H c s = — (Sl^ /A)|l)(0|+h.c, where 
fipt and A are the laser Rabi coupling and detuning, re- 
spectively. Therefore we see that an effective, small but 
non- vanishing, energy scale E\ := fij^/A is introduced. 
(Note that in order for perturbation theory to be valid 
one must have J7r <C A, which in turn implies E\ <C A.) 
Yet another energy scale is provided by the spectral width 
E2 of the laser pulse shape fiu,(i); in order to suppress 
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unwanted real transitions, one must impose in addition 
that E 2 <C A. At any rate, the appearance of these new 
system-energy scales implies that once again the SCL- 
type contradiction applies. On the other hand, we can 
make both E\ and E 2 small at the price of lengthening 
the gating time (r g ~ max{l/i?i, I/E2}). This implies, 
once again, an adiabatic limit and the applicability of 
the WCL. Therefore it appears that as long as one re- 
stricts manipulations to adiabatic ones (thus contradict- 
ing Al), quantum computing with degenerate qubits is 
possible even in the Markovian limit. We expand on this 
viewpoint below. 



G. Impure Ancillas 

With respect to the SCL: "Do ancillas really need to 
be pure?'''' 

What precisely is the role of the ancillas in QEC? A 
popular answer is that they serve as an "entropy sink" 
for the errors accumulated during the quantum compu- 
tation. This entropy in the system arises from the en- 
tanglement between system and bath, and the role of 
the ancillas is to remove this entanglement. I.e., in a 
perfect quantum error correction step the entanglement 
between system and bath is transferred to the ancillas 
and bath. A natural objection to our SCL-based incon- 
sistency is to claim that, in fact, ancillas need not be 
pure, or could perhaps even be highly mixed. However, 
this is not supported by the (current) standard theory of 
FT-QEC. Consider, e.g., an error correction circuit based 
on the Steane 7-qubit code. It takes as input ancillas pre- 
pared in the \tp) a = (|0l) + | 1_l) )/ a/2 state, where |0l) 
and |1l) are codewords. The physical qubits which com- 
prise such ancillas, are coupled bitwise via CNOT gates 
to the physical qubits making up the encoded data qubits 
in the circuit. If instead we input an ancilla in a mixed 
state, this is equivalent to inputting a classical mixture 
with erred codewords, e.g., (1 — p)\tp) a (ip\ + p\(j>) a ((j>\, 
where \4>)a is an erred codeword. If one of these errors 
is a phase-flip, it feeds back (via the CNOT gates) into 
the data qubits, producing an error [g^. Without fault- 
tolerance this means that there are now two errors (in the 
ancillas block and the data block), which may be more 
than the code can handle. In FT-QEC theory such er- 
rors are accounted for, but their magnitude is bounded 
from above (e.g., p in the above example must be small). 
We note that an ancilla which is initially entangled with 
the data qubits (violating the assumption of being intro- 
duced into the circuit in a tensor-product state) is es- 
sentially equivalent to the case of an impure ancilla just 
described (tracing over the data qubits yields an impure 
ancilla state). 

A more general approach showing the importance of 
the assumption of pure ancillas is the following (fairly 
standard account of QEC). 
i) Preparation- 

Let the initial state of system + reservoir + ancillas, 



with respective Hilbert spaces Hs,'Hr,'Ha, be: p% RA = 
\^s){^s\ <8> |0fl)(0fl| <g> p A , where we have allowed for an- 
cillas in a mixed state pa- 

ii) System-reservoir interaction (decoherence).- 
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where e's denote the errors belonging to the set £ that 
the code C can correct, and where \e R ) are the corre- 
sponding states of the reservoir. The error operators U e 
are assumed to be unitary and with linear span of dimen- 
sion \£\. 

iii) System- ancilla interaction (syndrome extraction) .- 

This interaction takes the form Us a — J2 e es He <8> T e 
where the T e 's are unitaries over TIa such that T e |0^) = 
\e A ) and U e S I c <g) |e)(e| 13 : 
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= Yl U e \i) S )^s\Ui,®\e R )(e' R \®T e p A T el . 

(54) 

iv) Error recovery. — 

Unitary recovery is implemented via Usa — 
l^l -1 ^ 2 J2ee£ U\®Ir® \&a){^a\, where for unitarity we 
need (eA\e' A ) = 5 e ,e'- By applying Usa and tracing over 
both R and A (assuming the |e^) 's too are orthonormal) 
one obtains 



out 
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(55) 

In the case of a pure ancillas p A — lO^O^I one has 
(fA\T e p A T^\f A ) = \(f A \e A )\ 2 = 6f, e and therefore the 
ideal case p A ut = |'^s)(V's| is recovered. One can also 
consider the fidelity 

F := tyslpT'm 

e,f£E 

(56) 

Provided the error operators Uf satisfy the condition for 
a non-degenerate code (ij) s\U^U e \tp s) = <>/,e [3l|], one 
obtains F = J2 e e£( eA \TePAT^\e A ) = (0 A \p A \0 A ). 

Clearly, F = 1 iff p A = |CU)(0a|, i-e., the ancillas are 
pure. One can also consider non-unitary recovery via an- 
cilla measurements and conditional unitaries, with Kraus 
operators given by A e = \E\~ l/2 Ul <g> I R <g) |e,t)(e,4|- The 
conclusion that the ancillas' state must be pure is un- 
changed. 



13 We know that Hs = C®S<SiT> [<S=syndrome subsystem, dimP : 
\£\; Z>=remainder (=0 for subspace- based codes)] I3lll66ll67l . 
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We note that FT is obtained by adding concatenation 
and, in steps iii) and iv), preparing and coupling encoded 
ancillas with the system in a suitable way, e.g., as in the 
Steane-code example above. In this case it is permissible 
to allow slightly impure ancillas, and relax the assump- 
tions that, in step ii) the environment couples only to 
the system, and in steps hi,iv), the environment does not 
act. This formulation, however, does not allow arbitrarily 
mixed-state ancillas, as argued in the Steane-code exam- 
ple. While such a formulation of FT-QEC theory might 
still emerge (for example, by using algorithmic cooling 
techniques (63, |(5!|, which, however, at present assume 
perfect gates), it does not appear possible at present to 
relax the assumption of cold ancillas. 



H. Hot QC, cold ancillas, and fast QC-ancilla 
interactions in the SCL 

With respect to the SCL: "One can keep the ancillas 
coupled to a separate cold bath and then couple them for 
only a short time to the QC: what matters then is the T\ 
timescale and that one can be very long compared to the 
required ancilla-QC coupling time". 

Let us paraphrase this objection. If one can make Hsa 
(system-ancillas) very large then one could beat the rate 
of ancilla heating by strongly coupling the QC and ancil- 
las. I.e., suppose one would like to bring the ancillas in 
from their cold reservoir to couple to the system, which 
is coupled to a hot reservoir as required for the SCL. The 
ancillas then heat up fast, but there is a timescale associ- 
ated with this heating ("Ti"), which one wishes to beat. 
Now if one could make the system-ancilla coupling very 
strong then one could, presumably, use the ancillas (e.g. 
for syndrome extraction) faster than their heating rate, 
while they are still sufficiently pure for fault tolerance 
purposes. 

The simplest argument against this objection is the 
following. In the setting of the objection, the QC is de- 
scribed by the SCL (high T) while the ancillas are de- 
scribed by the WCL (low T). Strong and fast coupling 
between the QC and the ancillas is unacceptable accord- 
ing to the WCL because it is fast (only adiabatic manip- 
ulations are allowed) , and according to the SCL because 
it is strong ("strong" refers to the system's Hamiltonian 
part, while in the SCL this Hamiltonian is weak in com- 
parison with the system-bath coupling) . 

However, one could go on to argue that the ancillas are 
a different species than the QC qubits, and in particu- 
lar have a different intrinsic (less dense) energy scale, so 
that they are at low T on the scale set by the QC qubits. 
In this case both ancillas and QC are described by the 
SCL. Then the problem with the objection is the follow- 
ing: recall that in the SCL (see Section fill A|) one must 
rescale Hsr and Har. as Hsr/c and Har/c respectively, 
where here R denotes the common reservoir the system 
and the ancillas are coupled to. The heating rate is pro- 
portional to the square of the coupling strength to the 



reservoir, i.e., to 1/e 2 , and hence diverges in the SCL. 
Therefore to beat the ancilla heating process via fast 
manipulation of the system-ancilla coupling one would 
have to rescale Hsa at least by 1/e 2 , but this contradicts 
the SCL derivation, where in fact one must keep Hsa 
fixed while rescaling Hsr. The reason for this is that, in 
the SCL derivation, it is the system (now including the 
ancillas) that sets the timescale against which reservoir 
correlations must be accelerated. 14 



VI. ALTERNATIVES TO MARKOVIAN FT-QEC 

A. Nature of the non-Markovian errors in the 
WCL 

While we have pointed out that, in the WCL, the 
application of fast gates is likely to violate the condi- 
tions required for Markovian dynamics to persist, we 
have not been specific about the type of non-Markovian 
effects that will emerge. It is well known that FT- 
QEC is capable of dealing with errors that change due 
to the application of gates. Namely, assume (slightly) 
faulty gates correcting a specific error model described 
by a CP map A [recall Eq. JJJ], are applied in sequence, 
AU' N AU' N _ 1 ■ ■ ■ AU[, and these gates are (in some appro- 
priate norm) close to the ideal gates {Z/i}^!.-^ i.e., for all 
i, \\UlUj — I\\ <C 1. Then by inserting lljUi's everywhere 
one obtains the new sequence AnUnAn-iUn-i ■ ■ • AiUi, 
where A, := AU-Jjj , and FT-QEC is capable of deal- 
ing with such a (gate-modified) error model. However, 



Let us also consider the issue from the perspective of thermody- 
namics. This is not really necessary, since the arguments above 
about the SCL are rigorous, but is interesting in its own right. 
First, we remark that error correction should really be made to 
work at the common lower (initial ancillas') temperature. Heat- 
ing a part of a QC only to be closer to the Markovian limit is 
a suboptimal strategy, because it increases the strength of the 
noise and stimulates entropy production. Second, in standard 
FT-QEC heat (entropy) flows from the QC to the ancillas only, 
while in reality one should expect a flow in both directions and 
additionally an entropy production. To see this let us ignore for 
the moment the coupling of the QC to the bath, and consider 
ancillas coupled to a heat bath at temperature T. The ancillas 
can be kept pure by maintaining an energy gap 2> kT. Assume 
that the initial state of QC (C) and ancillas (A) is a product 
state \ipc) ® I iPa)- Switching on the interaction He a we induce 
an equilibration process (because the dynamics is Markovian) 
of C + A towards the Gibbs state pcA = &cp(—HcA/kT)/Z, 
which is entangled (here for simplicity He A contains not only 
the interaction but is the total Hamiltonian of C + A). After a 
single step of error correction the total state of C + A can be 
modeled by (1 — p)U\ipc) ® \iPa)U^ + PPCA, where U is unitary 
and < p <C 1. Then we switch off the interaction with the an- 
cillas. Whatever we do next separately with C and A, we cannot 
eliminate the error due to the entanglement present in the term 
PPCA ■ This type of incorrectable error accumulates and destroys 
FT-QEC. This is the back flow of entropy from the ancillas bath 
to the QC, mentioned above. 



15 



the non-Markovian effects that arise due to the appli- 
cation of fast gates in the WCL, will in general not be 
describable by a simple time-local modification such as 
A -> AU(Uj. This can be worked out, e.g., using the 
methods of Ref. ■ 

In order to formulate consistent alternatives to stan- 
dard, Markovian FT-QEC theory, it seems useful to start 
with a Hamiltonian formulation. As the discussion below 
will illustrate, it appears that a hybrid approach will be 
necessary, which combines alternatives to standard QC 
with a new version of FT-QEC. 

B. Adiabatic Quantum Computing (AQC) 

We keep A2 and A3, discard Al (fast gates), and 
work in a purely adiabatic mode, thus permitting a con- 
sistent WCL. This may indeed be possible using the adi- 
abatic quantum computing (AQC) approach of Farhi et 
al. [52| . At present there is little understanding of 
the fault-tolerance of AQC. Some recent works explore 
AQC in the presence of decoherence and/or control er- 
rors [7(3, El El El, El El Indeed, the subject of the 
adiabatic approximation in open quantum systems has 
only very recently been addressed |4(j , and used to study 
AQC in open systems 74] . Error correcting codes for 
AQC were introduced very recently in EH> but the corre- 
sponding universal Hamiltonians involve many-body in- 
teractions (four and six-body for 1-local and 2-local er- 
rors, respectively). 

C. Holonomic Quantum Computing (HQC) 

Another possibility for keepin g A 2 and A3, and dis- 
carding Al, is provided by HQC [53|,|54|. HQC is an adi- 
abatic scheme which relies on Abelian or non-Abelian ge- 
ometric phases to implement quantum logic gates. Quan- 
tum information is encoded in a degenerate set of eigen- 
states of a Hamiltonian depending on a set of controllable 
parameters, e.g., external laser fields (recall our discus- 
sion of degenerate states above). When these are adia- 
batically driven along a suitable closed path, the initial 
quantum state is transformed by a non-trivial unitary 
transformation (holonomy) that is geometrical in nature. 
The key point is that the geometrical nature of the quan- 
tum holonomics is believed to render HQC inherently 
robust against certain kinds of errors. This alleged fault- 
tolerance has only recently been seriously begun to be ex- 
amined |77| : the emerging picture is that, while stability 
against decoherence must still be assessed, HQC seems 
to exhibit a strong robustness against stochastic errors 
in the control process generating the required adiabatic 
loops 78]. Moreover, in the adiabatic limit of Marko- 
vian dynamics it has been show that the geometric phase 
of a single qubit coupled to a magnetic field is robust 
against both dephasing and spontaneous emission (but 
not against bit flips) Nevertheless, since deviations 



from strict adiabaticity are inevitable, and adiabaticity is 
particularly challenging to satisfy in open quantum sys- 
tems 0, it is tempting to combine HQC with FT-QEC 
in order to address the performance of HQC in the pres- 
ence of decoherence errors. Alternatively, we note that 
a hybrid approach that seems to be rather promising is 
the embedding of HQC within a DFS 80] . This amounts 
to realizing a set of universal quantum gates, acting on 
a DFS, by means of non- abelian quantum holonomies. 
This strategy brings together the "best of two worlds" : 
the quantum decoherence avoidance virtues of DFSs and 
the fault-tolerance of the all-geometric holonomic control. 
It is possible that such an approach can be implemented 
for quantum information processing in, e.g., trapped ions 
and quantum dots. 



D. Topological Quantum Computing (TQC) 

A robust way of performing quantum computations 
is based on excitations with fractional statistics, since 
they have several fault-tolerant properties built in. This 
idea is known as topological quantum computing (TQC) 
[55l IBH l57| . Physical realizations of the simplest versions 
of TQC have been considered in the literature, using, 
e.g., rotating Bose-Einstein condensates jU and super- 
conducting circuits [82^. Let C denote the manifold of 
quantum codewords. Strikingly, in TQC, one can have a 
trivial Hamiltonian, e.g., H\c = 0, but nevertheless ob- 
tain non-trivial quantum evolution due to the existence 
of an underlying topological global structure (boundary 
conditions). Quantum encoding is typically performed 
in a properly designed degenerate ground state C. This 
fact implies, for low enough temperature, an exponential 
suppression of errors on encoded quantum information 
due to thermal fluctuations. More importantly, topologi- 
cal features can render such a ground state stable against 
errors represented by local operators, namely error opera- 
tors that do not involve a number of qubits of the order of 
the size of the system. For example, in the so-called toric 
codes [5^, |H, qubits are encoded in the ground-state 
manifold of a lattice of interacting spins in such a way 
that degenerate ground states are mutually connected 
only via high powers (scaling linearly with lattice size) of 
local operators. Thus, here the fault-tolerance properties 
are already built-in at the physical level. However, while 
one can argue that topological encoding provides a stable 
and passive quantum memory, it is not self-correcting as 
in today's "effectively naturally fault-tolerant" classical 
architectures (see Ref. [84| for an eloquent exposition of 
this point). Moreover, it is important to realize that as 
far as we know, in its present state TQC still requires ac- 
tive intervention, in the form of FT-QEC, when one tries 
to compute fault-tolerantly. Indeed, Preskill writes in 
Ref. [83[p-62], "It is therefore implicit that the temper- 
ature is small enough compared to the energy gap of the 
model that thermally excited anyons are too rare to cause 
trouble, that the anyons are kept far enough apart from 
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one another that uncontrolled exchange of charge can be 
neglected, and in general that errors in the topological 
quantum computation are unimportant. If the error rate 
is small but not completely negligible, then the standard 
theory of quantum fault tolerance can be invoked to boost 
the accuracy of the simulation as needed" . Ref. j2£j takes 
this approach and explicitly lists A2 and A3 as necessary 
requirements for fault-tolerant TQC. In contrast, Al is 
definitely not required in TQC: one performs computa- 
tions by adiabatically dragging quasiparticles around one 
another, and these operations must be slow relative to the 
gap between the ground state and the first excited state. 
The larger the gap the easier it is to satisfy this adiabatic- 
ity condition, so this requirements is compatible with the 
thermal suppression of errors mentioned above. In addi- 
tion, TQC requires time-dependent controls to read out 
the encoded data (Ref. [8(| shows that all measurements 
can be postponed until the readout of the final result of 
the computation). However, a fully Hamiltonian analysis 
of the fault-tolerance of such measurements is still lack- 
ing. Nevertheless, one could argue that the error rate 
in a topological quantum computer could be made ar- 
bitrarily small by increasing the system size and careful 
engineering, so that (similarly to today's self-correcting, 
fault-tolerant classical computers), one could ultimately 
perform TQC without any active intervention other than 
read-out of the encoded data. An interesting, recent de- 
velopment in this direction was reported in Ref. [84| . 
which suggests that certain three-dimensional quantum 
spin-lattices might be self-correcting. 



E. Non-Markovian Quantum Computing 

We keep Al and A2 but discard A3 (the Markov ap- 
proximation) at least in part. This appears to be a rea- 
sonable approach in many cases, since the Markov ap- 
proximation is clearly a highly idealized limit (though it 
does hold remarkably well in some optical systems and 
in liquid state NMR). Indeed, the degree of accuracy to 
which the Markov approximation must be satisfied has 
been quantified, e.g., by Steane in [plj : the probability 
of an uncorrectable (i.e., non-Markovian) error per gate 
must be < 10 -10 for a computation involving 10 9 gates 
(this probability must scale with the input size, as ex- 
plained in Section IV D|) . Alternative approaches to deal- 
ing with non-Markovian baths are therefore of interest. 
For example, the papers present an extension 

of FT-QEC theory to a non-Markovian setting. We offer 
in this context the following observations: 

1. An important ingredient carried over directly and 
without change from Markovian FT-QEC theory, is 
the crucial role of the fresh and nearly pure ancillas 
IB1 I20I l2lj . We believe that the detailed mechanism 
for introducing and discarding ancillas at specific times 
should be reconsidered within a fully Hamiltonian frame- 
work. 

2. As recognized and discussed in |19| . the important as- 



sumption of a small norm of the system-bath interaction 
Hamiltonian (e.g., Eq. (58) in Ref. |2£j) is not satisfied 
for some standard models of open systems. For exam- 
ple, a linear coupling to a bosonic heat bath involves un- 
bounded interaction Hamiltonians and a high-frequency 
cutoff. In general, the assumption of a small norm of 
the system-bath interaction Hamiltonian is much stricter 
than the WCL and is not satisfied for most standard 
models of reservoirs. 

Another approach to fault-tolerance in a non- 
Markovian setting is the recently developed time- 
concatenated dynamical decoupling method | 22| (see also 
|87| for a version of dynamical decoupling with bounded- 
strength controls). However, comment 2. above about 
the small norm of the system-bath interaction Hamilto- 
nian applies here as well. Therefore more general meth- 
ods are required to deal with the full scope of baths one 
can expect in quantum computing implementations. A 
promising possibility in this direction is to incorporate 
fault-tolerant dynamical decoupling in a feedback loop. 



VII. CONCLUSIONS 

We have listed a set of minimal assumptions made 
in the theory of fault-tolerant quantum error correction 
(FT-QEC): 1) fast gates (on the timescale set by the in- 
verse of the relevant Bohr or Rabi frequency), 2) a supply 
of fresh and nearly pure ancillas, 3) a Markovian bath. 

We have also reviewed the only two known rigorous 
general limits leading to Markovian dynamics: the sin- 
gular coupling limit (SCL), which involves taking a high 
temperature limit, and the weak coupling limit (WCL), 
which requires either a constant or an adiabatic system 
Hamiltonian, and averaging over long times in compar- 
ison with the inverse of the relevant Bohr frequency. 
These two limits allow one to replace the reservoir au- 
tocorrelation function by a Dirac delta, which leads to 
the Markovian limit. 

A close examination of the assumptions of FT-QEC 
has led us to conclude that assumption 3 can be sus- 
tained together with assumption 1 in the SCL, and to- 
gether with assumption 2 in the WCL. However, it is 
not possible to maintain all three assumptions in cither 
the SCL or the WCL. We therefore conclude that, at 
present, there exists an inconsistency in the formulation 
of the theory of FT-QEC for Markovian baths. We have 
also listed a number of alternatives to Markovian FT- 
QEC which, from the point of view adopted here, are 
free of inconsistencies. However, none of these alterna- 
tives is so comprehensive as to include the full range of 
errors one might expect in a full-scale implementation 
of quantum computing. In particular, recent results on 
fault tolerance in non-Markovian settings [T^ . I20l . l2il.l22l . 
while representing a significant step forward, make a cru- 
cial assumption about the smallness of the norm of the 
system-bath interaction Hamiltonian, which severely re- 
stricts the class of physical reservoirs. 
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